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For systems described by finite matrices, an affine form is developed for the maps that describe 
evolution of density matrices for a quantum system that interacts with another. This is established 
directly from the Heisenberg picture. It separates elements that depend only on the dynamics from 
those that depend on the state of the two systems. While the equivalent linear map is generally 
not completely positive, the homogeneous part of the affine maps is, and is shown to be composed 
of multiplication operations that come simply from the Hamiltonian for the larger system. The 
inhomogeneous part is shown to be zero if and only if the map does not increase the trace of the 
square of any density matrix. Properties are worked out in detail for two-qubit examples. 



I. INTRODUCTION 

From the beginning, our understanding of quantum 
mechanics has involved both the Heisenberg picture 0,0 
and the Schrodinger picture 0, and the relation be- 
tween them 0, E3|. Full understanding has been for a 
quantum system that is closed, which means there is no 
need to consider that it might interact with anything 
else. An open quantum system is a subsystem S of a 
larger system and interacts with the subsystem R that is 
the remainder, or rest of the larger system (which could 
be a reservoir). The evolution in S is considered to be 
the result of unitary Hamiltonian evolution in the larger 
system of S and R combined. The Heisenberg picture for 
S is clear. A matrix that represents a physical quantity 
for S (an observable) is changed by the unitary trans- 
formation that changes every matrix that represents a 
physical quantity for the larger system. The Schrodinger 
picture for S has not been fully described. 

The state of S is generally not a pure state, even when 
the state of the larger system of S and R combined is 
a pure state, so there is no Schrodinger wave function 
and no Schrodinger equation for S. There is a density 
matrix that describes the state of S. We can expect 
BEIlHIIlEEQElIi that evolution in 
the Schrodinger picture for S will be described by linear 
maps of matrices applied to the density matrix for S. 
Attention has been focused on the particular case where 
the initial state of the larger system is described by a 
density matrix that is a product of a density matrix for 
S and a density matrix for R. Then the evolution in 
the Schrodinger picture for S is described by linear maps 
that are completely positive. They have been studied 
extensively E3, E3, EE IS El S 

Recently we considered the general case, where S and 
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R may be entangled in the initial state, so that the den- 
sity matrix for the initial state of S and R combined is not 
a product. We worked out examples for entangled qubits 
in some detail, and for any system described by finite 
matrices we showed how the evolution in the Schrodinger 
picture for S can be described by linear maps [l4|. 

There are observations to be made of properties and 
structure in what has been found. They show us new 
features of the Schrodinger picture for open quantum 
systems. Our first observation |b| was that the maps 
generally are not completely positive, and apply in lim- 
ited domains 0, 0] . Another observation [l5| is that 
putting the maps in an affine form |23l |25| , with homo- 
geneous and inhomogeneous parts, can separate elements 
that depend only on the dynamics from those that de- 
pend on the state of entanglement. This gives a picture 
that is simpler in some respects and may be easier to use. 
We develop that picture here, show directly how it relates 
to the Heisenberg picture, and find some new properties 
of the affine form. Then we work out a new larger set of 
examples for two entangled qubits. 

We do not consider equations of motion. A simple 
equation of motion would require that when the map 
that describes evolution for a time t\ is followed by the 
map that describes evolution for a time £2, the result is 
the map that describes evolution for the time t\ -\-ti- 
The maps for open quantum systems generally do not 
have this semi-group property. To assume that they do 
is to make an approximation [H HHH |H H3 EH . We 
want to simply describe the Schrodinger picture before 
considering approximations to it. 



II. FRAMEWORK 

We consider two interacting quantum systems S and 
R, both described by finite matrices: NxN matrices for 
S and MxM for R. We use the matrices F^ , F „ and 
F^v described in our previous paper 0. The F^ for 
fi = 0, 1, . . . N 2 — 1 are iV 2 Hermitian matrices for S such 
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that Fqq is Is, the unit matrix for S, and 

Tr s [FpoFvo] = N5^. (2.1) 

This implies that the F^q are linearly independent, so 
every matrix for S is a linear combination of the F^q . For 
example, the F M o for /z = 1,2,... A'' 2 — ! could be obtained 
by normalizing standard generators |32j of SU(N). The 
-Fo„ for v = 0,1,... M 2 — 1 are M 2 Hcrmitian matrices 
for R such that _Fbo is 1_r, the unit matrix for R, and 



(2.2) 



Every matrix for R is a linear combination of the Fq v . 
We use notation that identifies F^ with F m0 ® Ifl and 
F „ with lg ® F 0l/ and let 



F u 



F^®F 0v . 



(2.3) 



Every matrix for the system of S and R combined is a 
linear combination of the F^ v . 

We follow common physics practice and write a prod- 
uct of operators for separate systems, for example a prod- 
uct of Pauli matrices £ and S for the two qubits consid- 
ered in Section lYlIII simply as S3, not £<g)S. Occasion- 
ally we insert a <g> for emphasis or clarity. 

The matrices F^o for p, — 1,2,... iV 2 — 1 and F 0iy for 
v = 1, 2, . . . M 2 — 1 are generalizations of Pauli matrices 
(and like the Pauli matrices they have zero trace). We use 
them to describe density matrices the way we use Pauli 
matrices to describe density matrices for qubits. If II is 
a density matrix for the system of S and R combined, 
then 



n : 



1 



JV -1 



iV -1 M -1 



NM 



(2.4 i 



a=l a=0 0=1 

and the density matrix p for the subsystem S is 

1 ( N *~ 1 
p = Ty r U = - 1 + Y ( F oo)F a0 



(2.5) 



a=l 



so that 

(F a p) = Tr[F Q/3 n] (2.6) 

and in particular 

(F a0 ) = Tr s [i^o Tr^n] = Tr s [F a0 p] . (2.7) 
If U is a unitary matrix, then 

N 2 -l M 2 -l 

U^F tlv U= J2 Y ( 2 - 8 ) 

a=0 /3=0 

with the i M „ ; a p elements of a real orthogonal matrix, so 
^pAj-afi 1S taP;fj,v Since U^tU and UtU^ are 1, 



^00 ; a/3 — ^0q ^0/3 ; taf3 



on 



(2.9) 



III. AFFINE MAPS OF DENSITY MATRICES 

Suppose that in the system of S and R combined the 
matrices C that represent physical quantities are changed 
to U'CU by a unitary operator U . This is the Heisenberg 
picture. The mean values are changed to 

(C) u = (U^CU) = Tr [U^CU IT] = Tr [C UUU^] . 

(3.1) 

The result is the same if the matrices C are left un- 
changed and the density matrix IT is changed to UHW . 
This is the Schrodinger picture. 

Let A be a matrix for the subsystem S. In the Heisen- 
berg picture it is changed to U^AU so its mean value is 
changed to 

(A) u = (U*AU) = Tr s [A Ttr [UUU^] ] . (3.2) 

The Schrodinger picture for the subsystem S is that the 
density matrix p for S is changed to 



where 



P 

HQ) 
K 



Ttr [UILU*] = L(p) + K 



(3.3) 



Ttr 



u n 



(3.4) 



The L is a completely positive linear map that applies to 
any matrix Q for the subsystem density matrix or not. 
It has the property that L(t) is 1. The map L depends 
on U but does not depend on the state of R or on the 
state of entanglement of the subsystems S and R. 

The K is the only part of p u that can depend on the 
state of R or on the correlations between S and R. With 
the same K, the Eq. (|3.3|1 defines a map that applies to 
different density matrices p representing different states 
of S. The state of S can be changed without changing 
K. That is evident from Eq. I|3.4|l . Since IT and p ® i 
are density matrices that give the same mean values for 
any matrix A for S, 



Tr [AU] = Trs [^TtrIT] = Tr s [Ap] = Tr 



their difference does not need to change when the state 
of S is changed. Explicitly, from Eqs. I|2.4|l and (|2.5|l we 

see that 



AT 2 -1 M -1 



IT 



M NM 



(3.6) 



a=0 f3=l 



which does not depend on the (F a o) which describe the 
state of S. 

When we define a map, we consider all the (F a p) to 
be independent. The (F a0 ) describe the state of S. The 
(Fap) for P not are considered to be parameters of the 
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map that describe the effect of the dynamics of the larger 
system of S and R combined that drives the evolution of 
S. Different (F a p) for (3 not specify different maps. 
Each map applies to different states of S described by 
different (F a o). For each map there is one NxN matrix 
K. We explained this with examples in our previous 
paper [14(. We also mentioned there that an alternative 
map can be used in the special case of a product state 
where (F a p) is (F a o)(Fop); we will not consider that here. 

In the Schrodinger picture K accounts for the parts of 
mean values (A) u that in the Hcisenberg picture come 
from matrices U'AU not being matrices for S. With- 
out K, a mean value (A) u calculated in the Schrodinger 
picture would be 



« = Tr s [AL{p)\ 



Tr 



U ] AU p® 



Tr s 

1 " 

M 



A Tr, 



Up® — U f 
M 



(3.7) 



which is obtained in the Heisenberg picture by replacing 
IT with p ® jj, which cuts off the part of U^AU that is 
not a matrix for S. The full mean value {A) u is obtained 
by adding 



Tr s [ART] = Tr 



u ] au ( n - 



p 09 -r-r 

M 



(3.8) 



This equation l|3.8(l follows directly from Eq. (|3.4J) . In 
particular, we have 



Trs [F^K] = Tr 



M 



(3.9) 



for p = 1, 2, . . . iV 2 — 1, so, because Trs-K" is zero, 
K = Tt s [Ft**] F no 

N 2 -l 



E Tr 



u^f u0 u (n- P ®4 



M 



F, 



/MO- 



(3.10) 



This is how we actually calculate K, as in the examples 
for two qubits described in Section lVlIII We do not need 
to calculate UHU' for the whole system of S and R com- 
bined. We just calculate WF^qU for the basis matrices 
F M o for S, take the mean values of the parts that extend 
outside the matrices for S, and get K from Eq. (|3.1U|I . 



IV. PURITY DECREASE 

A property that depends simply on the presence or 
absence of K is that 



Tr 



(P U ) 2 ]<K[P 2 ] 



(4.1) 



for all density matrices p if and only if K is zero. Here 
is a proof. If K is zero then from Eqs. (|2.5|) . (|3.7|l and 



JV -1 



Hp) = jj i + E <^o)o^o 



(4.2) 



where 



(Fuofl = Tr 



1 
M 



JV -1 



E t^Q-ao(F a a) 



because Tr 



F a p P ^ ~m 



(4.3) 

is zero if [3 is not zero and i M o ; oo 



is zero when p is not zero. Let 

Xap = {F a0 ) for a = l,2,...N 2 -l, = (4.4) 
and x Q/ 3 = for other a, j3, and let 



JV -1 M -1 

a=0 /3=0 



(4.5) 



Then (F M o}o is X^ and 



Af-l JV J -1M J -1 

E «^) 2 ^ E E (*U 

/.( = 1 ^i=0 i/=0 

JV 2 -1 A/ 2 -l 

= E E 

q=0 /3=0 
N 2 -l 

= E<^o) 2 



(4.6) 



which implies the inequality Ij4.1|l . 

Suppose K is not zero. Then jOJ fails for at least one 
density matrix p. Let 



K=J2 X n\n){7 



(4.7) 



n=l 



P=^ = ^EWH- (4-8) 



AT 



For this p we have 
p y = L(p) + /v =^+F = E (4 + A «) l")H ( 4 - 9 ) 

n=l 

and, since 

E A» = TvK = 0, (4.10) 



JV 



JV 



4 



N 



Tr 



71=1 

N 



-Els**- 



71=1 v 7 n=l n=l V 7 

Tr [p 2 ] . 



(4.11) 



For this proof we assume that the inequality (|4.1|) holds 
when p is 1/iV. A map generally is meant to apply only 
to a limited set of density matrices p, where it repre- 
sents the result of the unitary Hamiltonian dynamics in 
the larger system of S and R combined. The examples 
worked out in Section lYlIII show there are maps that are 
not meant to apply when p is TL/N. In such cases, the 
assumption that the inequality (|4.1|) holds when p is 1 /N 
is a mathematical statement that does not have a direct 
physical interpretation. 



V. MAP OPERATIONS 

The map L can be done with multiplication operations 
simply related to U. Let 



M -l 



U= Y G(")Fo„ 

i/=0 

with the G(v) matrices for S. Then 

M 2 -l 

L(Q) = Y G{y)QG(v)\ 



(5.1) 



(5.2) 



u=0 



M -1 



Y G{^G(v) = jjTm [tfu] = t s (5.3) 



and 



M -1 



L{t s ) = Y G{v)G{u)^ = Jj TT R i Uuf ] = a s- (5.4) 



v=0 



Altogether 



M -1 



p u = Y G{v) P G{^+K. 



(5.5) 



The matrices G(v) depend on U and depend on the choice 
of basis matrices Fq v . Making that choice to conform 
with U can simplify the set of matrices G(v), as the ex- 
amples described in Section lVlIII will show. The matrices 
G(v) do not depend on the state of S and R. They can 
be calculated from U and used for any states. 



VI. LINEAR MAPS OF MATRICES 

We fill out the Schrodinger picture with a linear map of 
matrices Q for S that gives p u when applied to a density 
matrix p. It is 



Q^Q' = L(Q) + KTtQ 



or, in terms of the basis matrices, 



%' = 1 + NK 



L(F a0 ) 



(6.1) 



(6.2) 



for a = 1, 2, . . . iV 2 — 1. This is the only linear map that 
can give p u for a variety of density matrices p described 
by Eq. 1)2.5(1 . Since K is the same for all p, it cannot 
come from the terms with variable coefficients (F a o). It 
can only be part of l'. 

We described this map in our previous paper 01 . We 
approached it differently there. We considered first the 
map of mean values (F a o) for the basis matrices for S and 
then the consequent maps of density matrices p and of 
the basis matrices 1 and F a Q. By working out examples 
of two entangled qubits, we found that this linear map 
(|6.1(l is generally not completely positive, that there is a 
limited domain in which it maps every positive matrix to 
a positive matrix, and that there is a limited domain in 
which it represents the effect of the dynamics of the larger 
system. We call these domains the positivity domain and 
the compatibility domain. We considered the description 
of the linear map l|6.1|l by 



N 



Qrs ^ F> r j- s kQjk 
3,k=l 



(6.3) 



and by 



N- 



Q' = Y C(n)QC(n? - Y G(n)QC{n)l (6.4) 



71=1 



Here we only make two comments. From Eq. (16.111 we 
see that K contributes K rs 8jk to B r j :s k in Eq. 1)6.3(1 . 
From Eqs l|6.1(l and (|5.2(l we have 



M —X 



Q'=Y G{v)QG{ V y +KTyQ. (6.5) 



In our experience with examples, the matrices G{y) and 
K in this equation have been significantly simpler than 
the matrices C(n) in Eq. (|6.4|) . 



VII. QUANTUM PROCESS TOMOGRAPHY 

How is such a map found? Is it observable? What can 
be seen in experiments? Is the map determined |23| by 
the effect of the dynamics on different density matrices 
pi It is if the compatibility domain contains an open set 
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of values for the (F a o). Then for each a from 1 to TV 2 — 1 
there are states of S, with density matrices p described 
by Eq. 1)2.5(1 . that differ only in the value of {F a0 ) for 
that one a. Between two of these states, the difference 



We let 



1 3 

K = 5E K i E i 



(8.3) 




so 



is just 



(F a o)Ko 



Ap u = -(A{F a0 ))F^ . 



(7.1) 



(7.2) 



This determines F' a0 . 
these F' aQ for a from 1 to N 2 



The map is specified by 1' and 
1. When all the F' a0 are 
known, 1' can be found from any p u . If the compatibility 
domain contains the state where (F a o) is zero for all a 
from 1 to TV 2 — 1, so p is •jkl, then l' is determined by 



p u = -1' 



(7.3) 



for that state, but as examples described in Section lYlIII 
show, the compatibility domain does not always include 
that state. 

The compatibility domain is the set of density matrices 
p that can be affected by the dynamics for the states be- 
ing considered for the larger system of S and R combined. 
If there are enough density matrices in the compatibility 
domain that are accessible to experiments, the map can 
be determined experimentally. The choice of the density 
matrices p to be used |2^| will depend on the particular 
situation. Density matrices that are handy for one sit- 
uation may not even be in the compatibility domain for 
another situation. 



VIII. TWO-QUBIT EXAMPLES 

We consider two qubits described by Pauli matrices 
Ei, E2, E3 for S and Hi, Ha, H3 for R, so Fjq is Ej and 
F k is Hfe, which implies Fjk is EjHfe, for j, k — 1,2,3. 
The density matrix for the two qubits is 

(3 3 
a + E(%)^ + ^(H fc >H, 
3=1 fc=l 



- E (EjHfc)EjSj. 

3,k=l 



(8.1) 



and the density matrix for S is 



Tr, 



= Tr s [EjK] 



(8.4) 



We write (E) for the vector with components (Ei), (E2), 
(E3) and k for the vector with components Hi, K%, K3, 
and write |(E)| and |k| for the lengths of these vectors. 
We write E for the vector whose components are the 
matrices Ei, E2, E3. 

Our examples are for different unitary matrices U. For 
the first set we think of U as describing the dynamics of 
the two qubits. In the second set U describes a Lorentz 
transformation of the spin of a massive particle for states 
with two possible values of the momentum. This illus- 
trates how the maps developed for dynamics can be used 
for other transformations as well. 



A. Interaction Hamiltonians 

The first examples are motivated by considering 
Hamiltonians that have only interaction terms, no free 
Hamiltonian terms, as in an interaction picture. From 
Ei,fe=i ijk^i-k we can get Ej=x7j E i s j b Y making a 
rotation in each qubit |33j| and redefining the Ej and H& , 
so to choose an example we let 



U 



_ e -i5(7iSiSi+7 2 S2H 2 +73S3=3) 



(8.5) 



(where 71, 72, 73 can be functions of time). The three 
matrices EiH 1; E 2 H 2 , E 3 H 3 commute with each other 
(The different Ej anticommute and the different Hj anti- 
commute, so the different EjH, commute). That allows 
us to easily compute 



-172 £ 2 



-273 £32:3 



= Ei cos 72 cos 73 + Hi sin 72 sin 73 

-E 2 H 3 cos 72 sin 73 + E3H2 sin 72 cos 73 

(8.6) 

using the algebra of Pauli matrices, and similarly 



U^^U = E2 cos 73 cos 71 + ^,2 sin 73 sin 71 

-E3H1 cos 73 sin 71 + E1H3 sin 73 cos 71, 



•7) 



J7^E 3 C7 = E3 cos 71 cos 72 + H3 sin 71 sin 72 

— EiH 2 cos 71 sin 72 + E 2 Hi sin 71 cos 72. 

(8.1 



6 



Interchanging U and U' has the same effect as changing 
the sign of every jj . Thus we see that 



= Ei COS72 COS73 
L(E 2 ) = E 2 cos 73 cos 71 
L(E 3 ) = E 3 cos7iCos7 2 . 



(8.9) 



Taking mean values in Eqs. ffTBjl . iPTTjl and (JE2I) gives 
the (Ej)' 7 , from which we see that 

Ki = (Si ) sin 72 sin 73 - (E 2 S 3 ) cos 72 sin 73 

+ (S 3 S 2 )sin7 2 cos7 3 



«2 



K.3 



(5 2 ) sin 73 sin 71 - (E3S1) cos 73 sin 71 

+ (EiS 3 )sin7 3 cos 71 

(5 3 ) sin 71 sin 72 - (EiS 2 )cos7i sin 72 

+ (E 2 Si) sin 71 cos 72. 



(8.10) 

We can construct K from Eq. (|8.3() and then get the 
linear map from Eq. ((6.1(1 or use Eq. 16.2(1 to get the 



.11) 



which determine the linear map. For the description of 
the linear map by Eq. ((6.3(1 we find that 



(1 + k 3 + G\Ci K\ — in2 

1 + Ks — C1C2 C2C3 — C3C1 ,. 1 , 

Kl + iK2 C 2 C 3 -C 3 Ci I-K3-C1C2 " ' 1 

C 2 C 3 + C 3 Ci ki + zk 2 

where = cos7 2 ; for i = 1,2,3, and the rows and columns of the matrix are in the order 11, 12, 21, 22. You can 
check that this is correct because it does give E 1; E 2 , E 3 , 1' that agree with Eq. ((8.111) . 





FIG. 1: (color online) The compatibility domain for (Hi), (S2}, (S3}, (E1H2}, (E1H3), (E2S1), (E2H3}, (E3S1} and (E3H2) all 
equal to \. (A) The whole compatibility domain inside the unit sphere. (B) The section of the domain in the (Ei, E2) plane. 
The sections of the domain in the (Ei, E3) and (E2, E3) planes are identical. 



The examples described in our previous paper [1J| are 
obtained as a particular case by letting 71 and 72 be zero, 
taking 7 3 to be u>t, and changing the S 3 in U to Ex. The 
new set of examples is much richer. In U there are three 
parameters instead of one. In K there are nine mean 



values: the three (Sfc) and the six (Ej-E/.) for j ^ k. The 
maps described in our previous paper depend only on 
(ExSx) and (E 2 Sx). 

Different values of the (Ej.) or (E^Sfc) in K generally 
give different maps. Each map is made to be used for a 
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particular set of states described by a particular set of 
density matrices p, or a particular set of (E), which we 
call the compatibility domain. It is the set of (E) that are 
compatible with the (Sfc) and (EjSfc) in K in describing 
a possible initial state for the two qubits. The increased 
number of (Ek) and (Ej-Sfc) in K means that the com- 
patibility domains are more restricted and varied. It is 
difficult to describe general features of the compatibility 
domains beyond the fact that they are convex [Tj|. 

In a larger domain, which we call the positivity do- 
main, the map takes every positive matrix to a positive 
matrix. The positivity domain is the set of (E) for which 
< 1. It depends on both the jj in U and the (S&) 
and (EjSfc) in K 7 so the variety of positivity domains is 
larger than the large variety of compatibility domains. 
We have looked at several examples. 

These examples exhibit new features. For the exam- 
ples described previously 0], the compatibility domain 
is not changed by reflection through the origin in the 
space of the (E); if (E) is in the compatibility domain, 
then so is —(E). The origin, the zero (E), is always in the 
compatibility domain. We can see from Figs. 2]El an dEl 
that these properties do not hold as a general rule. An- 
other property of the examples described previously Q 
is that the compatibility domain is the intersection of all 
the positivity domains for the same values of the 
and (EjSfc) in K. We can easily see that this also is not 
generally true. There are simple cases where the zero (E) 



is in every positivity domain but not in the compatibility 
domain. 

For example, suppose (Si) and (E3S1) are positive and 
all the other (S^) and (E^S^) for j ^ k are zero. Then 

ki = (Si) sin 72 sin 73, 

k-2 = -(E 3 Si)cos7 3 sin7i, 

k 3 = (8.13) 

and 

\k\ 2 < sin 2 73 + cos 2 73 = 1. (8.14) 

This implies that the zero (E) is in all the positivity do- 
mains for different 71, 72, 73, because (E) 17 is Ft when 
(E) is zero. 

We can see from Fig. [SJtliat for cases of this kind there 
are compatibility domains that do not contain the zero 
(E). We can easily show that there are many such cases. 
First we show that there is a substantial compatibility 
domain for any values of (Si) and (E3S1) short of the 
limit where (Si) or (E 3 Si) is 1. We find a set of (E) for 
which the matrix 

II = i(l + (E 1 )E 1 + (E 2 )E2 + (E 3 )E 3 

+ (E 3 Si)E 3 S 1 + (S 1 )Si) (8.15) 




FIG. 2: (color online) (A) Sections in the (£1 , E2) plane of the compatibility domain (dotted region) and the positivity domain 
(thick red curve) for (Si), (H 2 ), (H 3 ), (EiH 2 ), (E1E3), (E 2 Ei), (E 2 E 3 ), (E3E1) and (E 3 E 2 ) all equal to \ and 71 = 2\/5, 
72 = 2\/3, 73 = 2^2. (B) Sections of where the map takes the compatibility domain, the positivity domain and the unit sphere 
(thick red curve plus the dashed blue curve). The dotted circle is the section of the unit sphere 
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FIG. 3: (color online) The compatibility domain for (si) and (E3E1} equal to -7- and (H2), (H3), (E1H2), (E1H3), (E2S1), 
(E2H3}, (E3H2} all equal to zero. (A) The whole compatibility domain inside the unit sphere. (B) The section of the domain 
in the (Ei,E3) plane. (C) The section of the domain in the (E2,Es) plane. The compatibility domain does not intersect the 
(Ei, E2) plane. 



is positive (which implies that it is a density matrix and 
the set of (E) is in the compatibility domain). Let 

n= i(l + <5i)3i+M). (8.16) 

We see that M commutes with Si. There is a basis of 
eigenvectors of Si and M that diagonalizes EL From 

M 2 = (Si) 2 + (E 2 ) 2 + (E 3 ) 2 + (S 3 S!) 2 + 2<£ 3 )(£ 3 £i)£ 1 

(8.17) 

we see that the magnitude of the eigenvalues of M is 

V / (Si) 2 + (S 2 ) 2 + ((E3)±(S3S 1 }) 2 (8.18) 

when the eigenvalue of Si is ±1. The eigenvalues of II 
are all non-negative if 

(Si) 2 + (E 2 ) 2 + ((S 3 ) ± (S 3 Sx)) 2 < (1 ± (Si)) 2 . (8.19) 

These inequalities say that (E) is in the intersection of the 
two spheres of radii 1 ± (Si) with centers at =f(EsSi) on 
the 3 axis. There is a substantial intersection, so there 
is a substantial compatibility domain, for all values of 
(Si) and (E3S1) short of the limit where (Si) or (E3S1) 
is 1. For example when (Si) and (E3S1) are equal, the 
intersection is just the smaller sphere. 

Now we show that the zero (E) is not in the compat- 
ibility domain when (Si) + (E3S1) is larger than 1. We 
show that then there is no density matrix 

n = - (1 + (EiSi)EiHi + (E 2 S 2 )E 2 S 2 + (E 3 S3)E 3 S3 

+ (E 3 Si)E 3 Si + (Si)Si) (8.20) 
because then no matrix of this form can be positive. Let 

n = - (1 + (E 2 S 2 )E 2 S 2 + (E 3 S3)E 3 S3 + (Si)Si + W) 

(8.21) 



where 

W= [(EiS 1 )Ei + (E 3 Si)E 3 ]S 1 . (8.22) 

Because W 2 is (E1S1) 2 + (E3S1) 2 , the eigenvalues of W 
are 

± v / (EiSi) 2 + (E 3 S 1 ) 2 . (8.23) 

The matrices Si and W commute and make a complete 
set of commuting 4x4 matrices. Their eigenvalues label 
basis vectors for the four-dimensional space, one for each 
of the four combinations of the two eigenvalues of Si 
and the two eigenvalues of W. In particular, there is a 
nonzero vector ^ that is an eigenvector of Si and W for 
the negative eigenvalues of both. Since ('F, E 2 S 2 \I/) and 
(^, E 3 S3\I/) are zero, it gives 

(*, n*) = i||*|| 2 (l - (Ei) v /(E 1 Si) 2 + (E 3 S 1 ) 2 ) 

(8.24) 

which is negative if (Si) + (E3S1) is larger than 1. 



B. Lorentz transformations of spin 

These examples are abstracted from Lorentz transfor- 
mations of the spin of a particle with positive mass and 
spin i for two possible values of the momentum [34| . Let 

l7 = D 1 i(l + Si)+i? a i(l-Hi) (8.25) 

with Di and _D 2 the unitary rotation matrices made from 
E, so that 

D\£Di =i?i(E), D\tD 2 = i? 2 (E), (8.26) 
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for rotations R\ and i? 2 ; each i?(E) is simply the three- 
dimensional vector E rotated by R. In the application 
to Lorcntz transformations of spin [34|], E describes the 
spin of the particle, Si is the Pauli matrix for states 
with two different momentum values p 1 and P2 that is 
+1 when the momentum is pi and —1 when the mo- 
mentum is P2, and R\ and R2 are the Wigner rotations 
for the Lorentz transformation for pi and p 2 . Then U 
describes the Lorentz transformation of the spin in the 
system of two qubits where one qubit is the spin and the 
other is made from the two values of the momentum |34| . 
We would not have thought to consider an example this 
simple had it not come to us in an interesting application. 
From Eqs. and ijQjj) we get 



i?i(S)i(l + Sx)4 
\ + J R 2 (E) 



ifc(E)i(l-Hi) 



i?i(E) -i? 2 (E) 



.27) 



L(E) = -[i? r 1 (s) + i? 2 - 1 (E; 



i?i(£)-ii 2 (£) 



G(0) = ~(D 1 + D 2 ) , G(1) = ~(D 1 -D 2 ) 



(8.28) 



.29) 



(8.30) 



The mean values (E) are mapped to 

(X) u = (tfZU) = i(i?i(E) + i? 2 (E))+« (6 

and the density matrix p of Eq. (|8.2|l is mapped to 
1 



.31) 



P 



1 + (E) C 



\ U + ^(i?i(E) + i? 2 (E)} -E + k-E 



1 



k • E 



with 



L(p) + X 
l(l' + (E).E' 



£' = L(E) , i' = 1 + 2K. 



.32) 



(8.33) 



If (i?i(E)Si) = (i? 2 (E)^i) then K is zero and the map 
is completely positive. In fact the map is the same as it 
would be if (ESi) were zero. If (SSi) is zero then 



4d -s i: 



.34) 



In the application [34| , the mean value of the spin is the 
same for both momenta p! and p 2 . 

If (i?i(E)Si) ^ (i? 2 (E)Si) there are positive matrices 



1 



(1 



that are mapped to matrices 



1 



(1 



.35) 



(8.36) 



that are not positive. To see this it is sufficient to consider 
the case where R\ is the identity rotation and R 2 is an 
arbitrary rotation R. The general case can be recovered 
by taking R to be i?^~ 1 i? 2 and joining the same rotation 
i?i onto both the identity and R^ 1 R 2 to restore i?i and 
R 2 . The overall rotation R\ will just rotate a! and not 
change the non-positive character of the matrix 
Hence we consider 



1 



R(a)} 



i?(E) 



(8.37) 



Let a be along the axis of R so that R(a) is a. Then 
1 



a = a + 



E-.R(E) 



Choose the direction of a so that 



^ 1 
a ■ — 
2 



i?(E) 



Si) >0. 



(8.38) 



(8.39) 



Then |a'| > \a\. When |a| approaches 1, the matrix 1 
is not positive. 

For these maps, K depends only on (ESi), not on 
(ES 2 ), (ES 3 ) or (S). The compatibility domain is the 
set of (E) that are compatible with specified (ES^ in 
describing a possible state for the two qubits. It is very 
similar to the compatibility domain for the examples de- 
scribed in our previous paper. More precisely, K depends 
011 



i? 1 ((ES 1 ))-i? 2 ((ES 1 )) = R 1 (jESi) - i? 1 - 1 ^2((ES 1 ))^ 

(8.40) 

Let the 3 axis be along the axis of i? 1 " 1 i? 2 . Then (E3S1) 
is not changed by R^ 1 R 2 , so it drops out, leaving only 
(E1S1) and (E 2 Si) in K, and the compatibility domain is 
the set of (£) that are compatible with specified (E1S1) 
and (E 2 Si) in describing a possible state for the two 
qubits. This is exactly the compatibility domain for 
the examples described in our previous paper |l4j |. The 
equations and drawings that describe the compatibility 
domain there (Eqs^ (2.58), (2.64), (2.75), (2.77) and 
Figs. 2 and 3 in [l4|) apply here as well. 

The positivity domain is not the same as for the ex- 
amples described in our previous paper |l4j . It is the 
domain in which every positive matrix is mapped to a 
positive matrix, or the set of (E) for which 



1 r 

2 



fli((E)) +J R 2 «£» 



+ K 



<E)< 



< 1. 



.41) 
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C. The size of the inhomogeneous part 

How big can \k\ be? For the Lorentz-transformation 
examples described in Section IVHI Bl we find that the 
limit on \k\ is 1, but that \k\ can have any value short of 
that limit |34[- We do not know if values of |k| larger than 
1 are possible for the interaction-Hamiltonian examples 
described in Section lYlII Al we have not found any values 
larger than 1. Here is a proof that \k\ cannot be larger 

than which is about 1.62. 

From Eqs. (|8.4|l and (|3.4|) we have 



Kj = Tr s 



Tr 



u [ n - p ® 4 ) 



1 



WEjU [n-p®- 



.42) 



with n and p the density matrices of Eqs. (|8.1|) and 
(|8.2(1 . Since the trace of the product of matrices has the 
properties of an inner product for the real linear space of 
Hermitian matrices, it follows that 



( K3 ) 2 <Tr [(*7+E 3 C/) 



Tr 



1 N 



i)(l>> 2 +E<^> 5 



(8.43) 



From 

1 > Tr [n 2 ] 



fe=i 



we get 



M 2 <3-^(S,} 5 



(8.44) 



.45) 



This holds for any unitary matrix U, so it holds when 
U is replaced by UV with V a unitary rotation matrix 
made from the Pauli matrices WTijU that rotates U^TJJ 
to V^WtUV so that the only nonzero component of k 
for is k 3 and \k\ 2 for U is (k 3 ) 2 for [/P. Thus 

we conclude that 



<3-|(S)| s 



(8.46) 



We consider mean values (Ej), (Sfc), (SjHj.) that de- 
scribe a possible initial state for the two qubits. Then 
also 



Tr, 



( U^EU) - Tr, 



XL(p) 



< 1 



because (lI^U^ < 1 and, as in Eq. I|4.6|l 



.47) 



3 3 

^(TrsP^p)]) 2 = E(^o)o) 2 

<E(^o) 2 -E(^> 2 . 

(8.48) 



As a function of | (E) | , the bound l|8.46|l decreases and 
the bound (|8.47|) increases. The two bounds allow the 
largest when they meet. Then |(E)| is ~ 1 ^" v ^ and 
the largest |k| allowed is • ^ s the limit of large |k| 

is approached, the room for variation in |(E)| decreases, 
so room for the compatibility domain decreases, as our 
examples have shown. 
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